Incomplete devil's staircase in the magnetization curve of SrCu 2 (B03)2 
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We report on NMR and torque measurements on the frustrated quasi-two-dimensional spin-dimer 
system SrCu2(BC>3)2 in magnetic fields up to 34 T that reveal a sequence of magnetization plateaus 
at 1/8, 2/15, 1/6, and 1/4 of the saturation and two incommensurate phases below and above the 
1/6 plateau. The magnetic structures determined by NMR involve a stripe order of triplets in all 
plateaus, suggesting that the incommensurate phases originate from proliferation of domain walls. 
We propose that the magnetization process of SrCu2(B03)2 is best described as an incomplete 
devil's staircase. 
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Coupled-dimer spin systems with singlet ground states 
have been extensively studied because of a variety of 
magnetic-field-induced quantum phase transitions. Mag- 
netic fields strong enough to close the zero-field energy 
gap for the triplet excitations may simultaneously induce 
a longitudinal magnetization and a transverse antifer- 
romagnetic order, breaking the U(l) rotational symme- 
try around the field direction. This corresponds to a 
Bose-Einstein condensation (BEC) of triplets [l|. In the 
presence of frustration the magnetization curve in higher 
fields may show plateaus as a consequence of Wigner 
crystallization of triplet bosons |2J . 

In studying the magnetization plateaus, the layered 
compound SrCu2 (603)2 has played a prominent role 
since the discovery of a sequence of plateaus at 1/8, 1/4 
and 1/3 of the saturation value [3] . Its crystal structure 
consists of two-dimensional layers of orthogonal spin-1/2 
dimers in a geometry known as the Shastry-Sutherland 
lattice 0, [H[ , in which frustrated interactions drastically 
reduce the kinetic energy of triplets, leading to the theo- 
retical predictions forplateaus having a superstructure 
of localized triplets [fj, [fl]. A symmetry-breaking su- 
perstructure of triplets has been indeed observed for 
SrCu 2 (B0 3 )2 by Cu and B NMR at the 1/8 plateau [I,©]. 

More recent experiments on SrCu2(B03)2 have re- 
vealed a rich phase diagram below the 1/4 plateau, 
but with controversial results. B-NMR [l(| and torque 
llj measurements have revealed at least two additional 
phases between the 1/8 and 1/4 plateaus and excluded 
any plateau with a smaller fraction than 1/8. Other 
torque measurements, however, have been interpreted 
as evidence of a sequence of plateaus at 1/9, 1/8, 1/7, 
1/6, 1/5, 2/9 and 1/4 [Hj]. State-of-the-art analytic ap- 
proaches [3, [H| to determine the effective low-energy 
triplet Hamiltonian have confirmed the early findings of 
plateaus at 1/3 and 1/2 0,0] and have predicted addi- 
tional plateaus at 1/9, 2/15, 1/6 and 2/9 according to 



Ref. [13j], and at 1/8, 1/6 and 1/4 according to Ref. [14|. 
These methods, however, are valid only when the ratio of 
inter- to intra-dimer couplings is less than 0.5, whereas 
in SrCu2(B03)2 this ratio is probably around 0.65. 

To resolve the discrepancies and to determine the mag- 
netic structures of plateaus, we have carried out torque 
and n B-NMR measurements in static high magnetic 
fields up to 34 T generated by a 20 MW resistive mag- 
net at LNCMI Grenoble. We observed plateaus at 1/8, 
2/15, 1/6, and 1/4 of the saturation and two incommen- 
surate phases below and above the 1/6 plateau. For all 
plateaus, the magnetic structures determined by high- 
field NMR involve a stripe order of triplets with differ- 
ent commensurability, suggesting that the incommensu- 
rate phases originate from proliferation of domain walls. 
These findings establish the first observation of an in- 
complete devil's staircase in quantum antiferromagnets. 

To understand the 11 B-NMR spectra, we recall that 
each boron site generates three NMR lines at the fre- 
quencies, f k = j(H + H int ) + kvQ, (k = -1,0, or 1), 
where vq = 1.25 MHz is the quadrupole splitting for the 
external magnetic field H applied along the c-axis (TBI 
and 7 = 13.66 MHz/T is the nuclear gyromagnetic ra- 
tio. The internal field H m t is produced by nearby Cu 
spins, 



Hi nt = J2Ms l c 



(1) 



where Aj is the hyperfine coupling constant to the z-th 
Cu site that has the local magnetization (S l c ) along the 
c-axis. We show in Fig. Q] a part of the low-frequency 
satellite spectra (k = — 1) at T = 0.43 K covering the 
most negative range of internal field (Hi nt < —0.13 T) 
attributed to those boron sites whose nearest neighbor 
Cu carries a large magnetization. For this range of H nit , 
the low-frequency satellite spectra do not overlap with 
the central line (k = 0) or high-frequency satellite (k — 
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FIG. 1. (color online) A part of the 11 B-NMR spectra covering 
the most negative range of internal field obtained at T = 
0.43 K in various magnetic fields. The values of magnetic 
fields correspond to the position of the spectral base line on 
the vertical axis. The purple, red, blue, and green spectra 
belong to the 1/8, 2/15, 1/6 and 1/4 plateaus, respectively. 



1), and therefore they directly give the distribution of 
Hi n t with typical accuracy of 2 mT. 

The spectrum at 27.9 T (purple) belongs to the 1/8 
plateau and exhibits sharp peaks that do not move in the 
entire field range of the plateau, features that are char- 
acteristic of a commensurate superstructure. The peak 
positions agree with previous reports (sl-fioj . Among the 
spectra displayed in Fig. [TJ we clearly identify two other 
ranges of field, 28.7 - 29.2 T (red) and 31.5 - 32.2 T 
(blue), in which the spectra present the same features, 
suggesting the existence of two additional plateau phases 
between the 1/8 and 1/4 plateaus. Outside these field 
ranges, the peaks are rather broad and their positions 
change continuously with the external field. 

The existence of new plateaus is also supported by 
the magnetization curve shown in Fig. [5] obtained by 
torque measurement using a cantilever technique. The 
non-coplanar structure of the CUBO3 layers allows an in- 
tradimer Dzyaloshinski-Moriya (DM) interaction, which 
has been shown to produce a transverse magnetization 
perpendicular to the magnetic field 16]. The torque 
(t) acting on the cantilever then consists of two terms: 
t = aM x H + 6M • VH, the first one proportional to 
the transverse magnetization and the second one to the 
longitudinal magnetization. Their relative size depends 
on the precise location and orientation of the sample and 
of the field profile inside the magnet, which are difficult 
to know. In the gapped phase of SrCu2 (1303)2 below 
15 T, where the longitudinal magnetization is strictly 
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FIG. 2. (color online) Inset: The magnetic-field dependence 
of the torque divided by field at T — 60 mK for the sam- 
ple positioned at ±10 mm off the nominal field center. Main 
panel: The thick black line represents the longitudinal mag- 
netization(see the text) with the vertical scale appropriately 
adjusted. The magnetization values at 1/8, 2/15, 1/6, and 
1/4 of the saturation are shown by the dashed lines. The 
horizontal bars indicate the field range of the plateaus de- 
termined by NMR. The open circles show the magnetization 
determined from the Cu-NMR shift data from Ref. 01. 



zero, t/H shown in the inset of Fig. [2] varies linearly 
with H due to the transverse magnetization. To elimi- 
nate this contribution and isolate that of the longitudi- 
nal magnetization, we took a linear combination of two 
measurements of t/H taken at different sample positions 
shown in the inset of Fig. [2j choosing their relative coef- 
ficients so that the resulting curve stays zero below 15 T. 
The result is shown as a black line in the main panel of 
Fig. [5J It agrees very well with the magnetization deter- 
mined from the Cu-NMR shift data below 26 T reported 
in Ref. Q (open circles). 

The magnetization curve shows a series of plateaus. In 
addition to the 1/8 plateau and the approach to the 1 /4 
plateau just outside the available field range, two other 
plateaus can be clearly recognized: one is adjacent to the 
1 /8 plateau and the other is approximately half way up 
to the 1/4 plateau. These field ranges agree perfectly 
with what we proposed from the field variation of the 
NMR spectra (the horizontal bars in Fig. [2]) . The mag- 
netizations of the first three plateaus scale as 1/8: 2/15 
: 1/6, which is partially consistent with the theoretical 
predictions in Refs. [13| and 1141 but does not agree with 
the values reported in Ref. [12j . 

Having established the sequence of plateaus, we now 
discuss the spin structure. The distribution of H lnt 
was obtained by an iterative method to deconvolute the 
quadrupole structure from the NMR spectra 17| . The re- 
sulting spectra are displayed in Fig. EJA) for all plateaus 
and two intermediate phases above and below the 1/6 
plateau. Let us recall that the spin density of one triplet 
is expected to be distributed primarily over three dimers 
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FIG. 3. (color online) (A) Deconvoluted 11 B-NMR spectra representing the distribution of Hi nt obtained at T = 430 mK unless 
explicitly indicated. The data for the 1/8 plateau at 27.5 T was taken from Ref. [9(. (B) An extended triplet and the nearest 
neighbor B sites. (C) The crystal structure of SrCu2(BOs)2. The arrows indicate the hyperfine coupling of a B site to the 
nearest and the second nearest Cu on the same layer (Al and A2) as well as the coupling to the third and fourth neighbors on 
the adjacent layers (A3 and A4). (D) The spin superstructure for the plateau phases. The thin black lines show the lattice of 
orthogonal Cu-dimers in one layer. The thick red lines show the triplet dimers carrying the largest magnetization mo in the 
same layer while the blue lines indicate these triplets on the neighboring layers. The unit cell of each superstructure is shown 
by the shaded area. 



as shown in Fig. [H, [l8| . The central dimer with a 

large magnetization (S c ) = too = 0.3 ~ 0.4 is surrounded 
by two orthogonal dimers with a negative magnetization 
m 2 = —0.1 ~ —0.2 (antiparallel to the external field) 
next to the central dimer and a positive magnetization 
mi = 0.2 ~ 0.3 at the ends. We call this cluster an ex- 
tended triplet. It is stabilized by the cooperative (non 
frustrated) action of the polarized spins of the central 
triplet on the two neighboring dimers, and by its strongly 
reduced coupling to surrounding spins due to frustration. 
Other dimers have much smaller magnetization. We de- 
fine Bo, Bi, and B2 as the B sites nearest to too, mi, 
and TO2, respectively. All plateau spectra show a group 
of lines near H int = —0.2 T isolated from other lines (the 
part shown in Fig. [I}, which comes from B and Bi sites. 
The magnetic structures of plateaus were determined as 
follows: 

(1) The coupling constant Ai is the sum of the short- 
ranged transferred hyperfine coupling and the classical 
dipolar coupling. The coupling to the nearest neighbor 



Cu spin (see Fig.[3fC)) is by far the largest. The second 
neighbor gives a very small contribution due to the can- 
cellation between the transferred and dipolar terms, and 
it is much smaller than the dipolar coupling to the third 
and fourth neighbors. Therefore, to zeroth order, H lnt is 
simply proportional to the magnetization on the nearest 
neighbor Cu site. 

(2) The zeroth-order NMR lines then shift and may 
be split by the dipolar field from the third and fourth 
neighbors, which are on the adjacent layers as shown in 
Fig. OUC). This means that the spectral features, which 
distinguish one plateau from the others, are primarily 
due to the interlayer stacking of extended triplets rather 
than the in-plane spin configuration. We found that the 
partial spectrum of Bo and Bi sites is specific enough 
in all cases to select a unique in-plane structure allowing 
at least one stacking pattern qualitatively consistent with 
the observed spectrum. Figure[3jD) shows the structures 
thus determined, including the one for the 1/3 plateau 
which we found by applying the same analysis to the 
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FIG. 4. (color online) Distribution of local magnetization in 
the 1/8, 1/6, 1/4 and 1/3 plateau phases deduced from the 
peak positions of the B-NMR spectra. 



spectrum reported by Stern et al. [19|. The unit cell of 
all these structures involve two layers as in the crystal 
structure. 

(3) As a consistency check, we have tried to account 
for the entire spectrum of the 1/8, 1/6, 1/4, and 1/3 
plateaus assuming the structure shown in Fig. [3jD). For 
these structures, we were able to assign a particular peak 
value of H int to each B site almost uniquely. The values of 
the local magnetization (5*) have then been determined 
by solving Eq. (1), since there is an equal number of in- 
equivalent Cu and B sites in a unit cell of the superlattice 
(see the Supplemental Material [20]). The resulting val- 
ues of mo, mi, and m%, displayed in Fig. HI change very 
little from one plateau to the next, as expected. 

All plateaus show stripe order of triplets as shown in 
Fig.|3][D). The structures for the 1/3 plateau and the 1/4 
plateau agree with the theoretical prediction @, The 
former structure is a simple close packing of extended 
triplets, while in the other structures they are progres- 
sively more spaced by sites having much smaller polariza- 
tion. The structure of the 1/6 plateau can be obtained 
from that of the 1/3 plateau by removing every other 
triplet from each stripe. Similarly, the square unit cell 
of the 1/8 plateau is obtained by removing every other 
triplet from each stripe of the 1/4 plateau. Note that 
this structure of the 1/8 plateau is different from the one 
with a rhomboid unit cell, which was proposed based on 
the previous Cu-NMR experiments [H (see Supplemental 
Material [13]). Finally, the structure of the 2/15 plateau 
consists of alternating domains of 1/8 and 1/6 structures. 
More precisely, three stripes of the 1/8 rhomboid struc- 
ture are followed by a single stripe of the 1/6 structure. 
This is a good example showing how the proliferation of 
domain walls leads to a structure of higher order com- 
mensurability. We should remark that the structures of 
the 1/6 and 2/15 plateaus are not the ones predicted 



theoretically [13|, presumably because the ratio of inter- 
to intradimer coupling is significantly larger than 0.5. 

Let us now discuss the spectra in the intermediate 
phases below and above the 1/6 plateau. As shown in 
Fig. |3][A), they share a common feature to the spectra 
in the plateaus, that is the Bo and Bi lines at large neg- 
ative values of H- lnt ~ —0.2 T well separated from the 
rest of the spectra. This indicates that the fully polar- 
ized extended triplets remain immobile at least within 
the time scale of NMR measurement (around 100 /is), 
i.e. they are completely localized. However, unlike in the 
plateaus, the spectra in the intermediate phases consist 
of a few broad continuous lines, which cannot be decom- 
posed into discrete narrow lines. Such a spectral feature 
indicates an incommensurate magnetic structure. This is 
quite natural since the density of triplets, i.e. the mag- 
netization, changes continuously, and therefore generally 
takes irrational values. Intermediate phases could have a 
commensurate structure if the system separates into two 
commensurate phases with their volume ratio changing 
with field, or if a commensurate fraction of triplets form- 
ing a superstructure coexists with the rest of delocalized 
triplets. However, the NMR spectra in such cases should 
have sharp peaks similar to the spectra in the plateau 
phases, clearly inconsistent with the experimental obser- 
vation. 

As an attempt to put the results in a broader perspec- 
tive, we conjecture that the sequence of phases revealed 
by the present measurements is the first example of an 
incomplete devil's staircase in the context of the magneti- 
zation curve of a quantum antiferromagnet. The concept 
was introduced in the investigation of commensurate- 
incommensurate transitions: a devil's staircase is an infi- 
nite sequence of commensurate phases with increasingly 
large commensurability built by the proliferation of do- 
main walls between domains of different commensura- 
bility [21]. It becomes incomplete when the infinite se- 
quences of high commensurability phases with small steps 
are replaced by incommensurate phases 

This interpretation is supported not only by 
the sequence of fractional magnetization ratios in 
SrCu2(B03)2, which is typical of a devil's staircase [2l| . 
but also by the stripe structure of the plateaus deter- 
mined by the present NMR experiments, which naturally 
allows the formation of domain walls, as seems to be the 
case in the 2/15 plateau. It would be very interesting to 
check this interpretation by a direct measurement of the 
wave- vector of the structure, which has become accessible 
for the required magnetic field range by neutron scatter- 
ing performed in pulsed field [23j. Finally, we note that 
it remains a challenge for theory to explain why higher 
commensurability plateaus (or even lower ones such as 
1/7) are unstable towards the formation of incommensu- 
rate phases. 
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